Incomplete pure dephasing of A^-qubit entangled W states 
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We consider qubits in a linear arrangement coupled to a bosonic field which acts as a quantum heat 
bath and causes decoherence. By taking the spatial separation of the qubits explicitly into account, 
the reduced qubit dynamics acquires an additional non-Markovian element. We investigate the time 
evolution of an entangled many-qubit W state, which for vanishing qubit separation remains robust 
under pure dephasing. For finite separation, by contrast, the dynamics is no longer decoherence- 
free. On the other hand, spatial noise correlations may prevent a complete dephasing. While 
a standard Bloch-Redfield master equation fails to describe this behavior even qualitatively, we 
propose instead a widely applicable causal master equation. Here we employ it to identify and 
characterize decoherence-poor subspaces. Consequences for quantum error correction are discussed. 

PACS numbers: 03.65.Yz, 63.22.+m, 03.67.Pp, 03.67.-a 
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I. INTRODUCTION 

In recent years, we witnessed great progress in the 
field of solid-state quantum information processing, like 
the coherent control of single qubitsii^i^ and two-qubit 
gatesi^ One of the major remaining challenges is decoher- 
ence: the interaction of the qubits with their environment 
reduces the indispensable quantum coherence and entan- 
glement of the quantum states. This relates to the scal- 
ability of the present few-qubit setups, because decoher- 
ence becomes more pronounced as the number of qubits 
increases. Understanding the scaling of multi-qubit deco- 
herence is also experimentally relevant, as many groups 
take the challenge of implementing more complex qubit 
architectures with solid-state devices. 

Not all many-qubit states are equally sensitive to the 
influence of an environment. Depending on the symme- 
tries of the qubits-environment coupling, there can exist 
distinguished subspaces of a 7V-qubit Hilbert space that 
are effectively decoupled from the environment and, thus, 
form so-called decoherence- free subspaces (DFS)i^i^iii^ 
These allow one to implement logical decoherence-free 
qubits with two or several physical qubits. Therefore, 
given an architecture for N physical qubits, it is essen- 
tial to identify the sets of robust quantum states that 
suffer least from decoherence. 

The substrate that supports the qubits also possesses 
its own degrees of freedom like nuclear spins and phonons. 
They generally are coupled to the qubits and, thereby, 
cause quantum dissipation and decoherencci^ d^i^^d^d'^d^ 
In experiments, one usually observes both dephasing and 
relaxation, with dephasing happening on a faster time 
scalei^'^^'^'' It has been noted that the reduced qubit dy- 
namics in principle can be solved exactly if the qubits ex- 
perience pure phase noise )^ d^d9i20i2ii22, 23^24^25 g^j^j ]^gj.g 

we focus on this case as well. In Ref. [2^, we recently 
presented explicit expressions for the dephasing of two 
initially entangled qubits. A central conclusion of that 
work is that the entanglement of a robust entangled 
state^^i^ is not perfectly stable but undergoes an ini- 
tial decay, stemming from the spatial qubit separation 
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FIG. 1: (Color online) Schematic representation of A'^ qubits 
in a linear arrangement. The qubits (green boxes) interact 
via a coupling to the substrate phonon field (red line) at their 
positions x^. 



sketched in Fig. [TJ This puts limitations on the applica- 
bility of the concept of decoherence-free subspaces, since 
at best decoherence-poor subspaces emerge instead. 

Here we discuss the effects of phase noise on the N- 
qubit generalization of the mentioned robust entangled 
state, namely the so-called W state;^ which is a coher- 
ent superposition of all states with exactly one qubit in 
state 1 while all the others are in state [cf. Eq. P^ ]. 
The W states play an important role in several pro- 
tocols for quantum information processing, for exam- 
ple quantum teleportation,— superdense coding^ and 
quantum gamesi^ In quantum optics, the W states 
have already been realized, first with three j^^'^^i'^^ and 
recently even with eight qubits."^" Proposals exist to 
produce W states in atomic gases'^® and in solid-state 
environmentSi^2i^i22 

In the limit of weak system-bath coupling, a successful 
and common approach to quantum dissipation is pro- 
vided by the Bloch-Redfield master equation.*" A cor- 
nerstone of this formulation is neglecting memory effects 
of the bath, so that one eventually obtains a Markovian 
master equation. The indirect interaction of two sep- 
arated qubits via the environment, however, introduces 
memory effects that arise when bath distortions can prop- 
agate from one qubit to another during a finite time. This 
timescale can be much larger than the intrinsic memory 
time of the bath. As will be detailed below, a direct 
application of the Bloch-Redfield approach to spatially 
separated qubits predicts spurious decoherence-free sub- 
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spaces. 

In the present work, we pursue a twofold goal: First, 
we consider an arbitrary number of N initially entangled 
qubits for which we present explicit expressions for the 
coherence loss. This shows how in a linear arrangement 
decoherence scales as a function of the number of qubits, 
and demonstrates the consequences of spatial qubit sep- 
arations. The stability of the W states with respect to 
the system size N has been studied previouslyj^ii^i^ 
however for local decoherence models where the qubits 
couple to effectively independent heat baths. Second, by 
making a weak-coupling approximation, we derive a non- 
Markovian master equation approach that captures the 
main effects of the spatial separation. This approach has 
the advantage of being more intuitive, while still allow- 
ing algebraic mcthods^'^ to be applied to the problem of 
decoherence. A comparison of the master-equation dy- 
namics and the exact dynamics therefore enables a better 
interpretation of the latter and a critical examination of 
the validity of the former. Moreover, the master-equation 
approach will be applicable as well to other problems that 
do not possess an exact solution. 

The paper is organized as follows: In Sec. [Ill we present 
our model, a system of N qubits in a linear arrangement. 
The qubits interact with a thermal bosonic field which 
causes decoherence of the iV-qubit state. In Sec. IIIIl we 
present for pure dephasing the exact time evolution of the 
concurrence for two qubits initially prepared in a robust 
entangled Bell-state and then generalize our results to W 
states of an arbitrary number of qubits. Starting from the 
general model, we derive in Sec.|TV]a master equation for 
the reduced dynamics, taking special care of the spatial 
qubit separation. Its solution is then compared with the 
exact results for pure dephasing. The derivations of the 
exact solution and of a convolutionless master equation 
are deferred to two appendices. 

II. QUBITS COUPLED TO A BOSONIC FIELD 

As a model for the N qubits in the bosonic environ- 
ment sketched in Fig. [1] we employ the Hamiltonian 

H = Ha+H^b, (1) 

where the qubits and the bosonic field in the absence of 
the coupling are described by the Hamiltonian^ '^°i^^ 

h ^ 

iJo = 2 X! ^'^^'^2 + X! fujjkb\bk ■ (2) 

v=l k 

The first term in Eq. ([2]) represents N qubits v = 
1, 2, . . . , with energy splittings Ml^, and Pauli matri- 
ces ■ Since we will not address the coherent control of 
individual qubits explicitly, the specific choice for the en- 
ergy splittings is not of major relevance. Note that there 
is no direct interaction between the qubits. The second 
term in Eq. ^ describes a bosonic field that consists of 



modes k with energies tuvk and the usual bosonic anni- 
hilation and creation operators hk and b\. We restrict 
ourselves to a linear dispersion relation Wfc = c\k\ with c 
being the sound velocity. 

Qubit V is located at position Xi, and couples linearly 
via the operator to the field, so that the coupling 
Hamiltonian reads 

N 

iJqb^ft^A,^,, (3) 

l/=l 

with 

?.=e(x.)-5]5fce^'=-''(6fc + 6L,) (4) 

k 

the bosonic field operator at position x^. We assume the 
coupling strengths to be real, isotropic, and identical for 
all qubits, i.e. gku = gk and g-k = gk- 

We choose an initial condition of the Feynman- Vernon 
type, i.e. at time t — to, the bath is at equilibrium and 
is not correlated with the qubits. Thus the total initial 
density matrix R{to) is a direct product of a qubit and 
bath density operator, 

Rito) = pito) ® ■ (5) 
Here, p is the reduced density matrix of the qubits and 




is the canonical ensemble of the bosons at temperature 
T and Z is the corresponding partition function. 

The dynamics of the qubits plus the environment is 
governed by the Liouvillc-von Neumann equation 

ih^^R{t)^[H^^{t),R{t)]. (7) 

The tilde denotes the interaction-picture representation 
with respect to Ho, i.e. A(t) = UQ{t)AUo{t), where 
Uo{t) = cxp{—iHo{t — to)/h}. We are exclusively inter- 
ested in the state of the qubits, so our goal is to find the 
time evolution of the reduced density operator p{t) = 
trbi?(i), where trb denotes the trace over the environ- 
mental degrees of freedom. In what follows, we will con- 
sider the density matrix elements Pm.n = the 
basis |n) = ^2, . . . , mat), where a^z\n) = (-l)"-|n) 
and Ui, = 0, 1. 

III. EXACT REDUCED DYNAMICS 

Generally, the coupling of a qubit to an environ- 
ment induces spin flips and also randomizes the rel- 
ative phase between the eigenstates of the qubit. If 
the coupling operator Hq^ commutes with the qubit 
Hamiltonian, the qubits experience the so-called pure 
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phase noise. Consequently, one finds [ffqb, Uq] = 
so that the interaction-picture qubit operators remain 
time-independent, X^{t) — X^. In particular, the time- 
evolution of the reduced density operator p is indepen- 
dent of the coherent oscillation frequencies Q.^, of the 
qubit s. 

In the following, we consider the coupling operators 
Xu = a^z which constitute a case of pure phase noise. 
The reduced qubit dynamics can then be solved analyt- 
ically. We defer the explicit derivation to Appendix \^ 
where we obtain 



with the amplitude dampingM 

/■°° , T, N 1 - cos(wi) , f huj 

Am,n(0 = / dujJ{uo) — ^coth 



(8) 



2fcBT 



N 



(9) 



and the time-dependent phase shift 
,. -, , — sin(wt) 

N 

X - (-l)"''+"''']cos(tjW)- 

v^u' — l 

(10) 

For ease of notation, we have set the initial time — 
and introduced the bath spectral density J(w) = 
— ck). The transit time of a wave between 
the qubits v and v' reads t^i,i — Xup'/c, where x^^i — 
\xu — ^u'l- We will focus on a linear arrangement of 
the N qubits and consider equal nearest-neighbor sep- 
arations Xy.jy+i — Xi2. Then the transit time becomes 

To elaborate on the impact of spatially correlated 
noise, we assume the chain of qubits to be embedded 
in a medium with a channel structure, i.e. we treat the 
bosonic field as effectively one-dimensional. Quasi one- 
dimensional geometries may, for example, be realized by 
carbon nanotubes or linear ion traps. They represent 
configurations in which the requirements of qubit pro- 
tection and adressability are well-balanced, and where 
we expect the effects of spatial noise correlations to be 
most overt. In this case, the spectral density is of the 
ohmic type,^ 



(11) 



wherein lo^ denotes a cutoff frequency which for a phonon 
field is the Debye frequency. For a more detailed discus- 
sion on the relation of the spectral density to the sub- 
strate geometry, in particular in the context of semicon- 
ductor quantum dots, we refer the reader to Refs. 25.45. 

For a physical realization with a GaAs substrate, the 
Debye frequency is of the order Wc = 5 x 10^^ Hz while 



the sound velocity is c = 3 x lO^ms"^. Then a qubit 
separation x\i = 100 nm corresponds to the transit time 
tyi — lO^'/ti'c. Likewise, for a temperature T — 10 mK 
we have k^T/TwOc — 10~^. This implies that lo'^^ is typ- 
ically the smallest timescale of the problem, while the 
transit time and the thermal coherence time can be of 
the same order. The present work is devoted to the ef- 
fects of spatial separation between qubits, but not of the 
finite extent of the qubits themselves. Note that in semi- 
conductor quantum dots, the finite width a of an electron 
wavefunction in the dot (typically below 100 nanometers) 
may lead to an effective cutoff frequency Wc ~ c/a that 
is smaller than the Debye frequency, but still larger than 
all other frequencies in the system4^i^i^ 

Both the damping ^ and the phases (|10p vanish at 
time t = 0, so that Eq. ^ is consistent with the ini- 
tial condition. As expected for pure dephasing, popu- 
lations are preserved, i.e. the diagonal matrix elements 
obey Pm,m(i) — /5m,m(0). This implies that generally 
neither the qubits nor the total system will reach ther- 
mal equilibrium. However, the relative phases between 
eigenstates will be randomized so that off-diagonal den- 
sity matrix elements — the so-called coherences — may 
decay, which reflects the process of decoherence. 



A. Dephasing of robust entangled states 

A most relevant decoherence effect in a quantum com- 
puter is the loss of entanglement between different qubits. 
In order to exemplify the impact of a spatial qubit sep- 
aration on decoherence, we consider as the initial state 
the robust entangled A^-qubit W state 
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|100...0)-l-|010...0)-l-...-l-|000...1) 



(12) 

For two qubits (A^ — 2) it has been shown that the two- 
qubit entanglement inherent in the state \W2), which is 
the symmetric Bell state, is robust under dephasing for 
vanishing spatial separation^^i^i while it decays for finite 
separation. ^-^^ 

Our motivation to focus on the initial states is 
twofold: First, W states play an important role in several 
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protocols for quantum information processing,- 
that their sensitivity to an environment is relevant in it- 
self. Second, among all fully entangled A^-qubit states the 
W states are special in that they maintain their A^-qubit 
entanglement under collective dephasing (i.e. for vanish- 
ing qubit separations). Now if already the W states start 
to loose their entanglement due to a finite spatial separa- 
tion of the qubits, then this is a strong indication that for 
other fully entangled A^-qubit states, the situation would 
be worse. Or, to put it simply, we wish to give the most 
optimistic estimate about A^-qubit decoherence and to 
our knowledge the best way to do that is by focusing on 
the W states. 

The fact that no bit fiips occur under pure dephasing 
is refiected in the structure of the exact solution (HI) : All 



4 



density matrix elements that are initially zero remain 
zero, so that for the state IW^at), the dissipative quantum 
dynamics is restricted to the states 

|j) = |00. .0), J = 1,2,. ..,7V. (13) 

Thus at most N"^ out of 2^^ density matrix elements are 
nonvanishing. Initially they are all equal, i.e. pjji (0) = 
1/N. From the Hamiltonian ^ it directly follows, 
that the states |j) possess the eigenenergies hwj — 

■| X^^^i ~ ^^"i ^ back-transformation of the co- 
herence pjji (t) to the Schrodingcr picture provides the 
phase factor exp[j(u;j — ujj')t]. 

a. Frequency shifts. One effect of the coupling to a 
heat bath is a frequency shift 5VLj which we obtain in the 
following way: Upon noticing that one can separate the 
phases (jlOl) into terms that depend on only j or j' , we 
write ipjj'i^) = 4'j{t) ~ 4'j'{t)- Each (j)j{t) turns out to 
consist of a finite contribution and a contribution that 
grows linearly in time, i.e. (f>j(t) = (pj(t) — dfljt. The 
latter leads to a (static) frequency shift. For the ohmic 
spectral density PT|) . we obtain 

N 

U^f' — l 

N 

"Pjit) = -f (-l)''^"^''^"' arctan[wc(t ± U^')] ■ 

(15) 



Thus the effective energy splitting of qubit j becomes 
h{VLj+5V!.j). Note that both i5rij- and(^j(t) depend on the 
transit times t^^i and the system size A'' but not on the 
temperature. They can be interpreted as a result from an 
effective coherent interaction of the qubits mediated by 
the vacuum fluctuations of the bosonic fieldj^^ where SVlj 
arises from an induced static exchange interaction and its 
onset is described by (pj{t). Note that the dominant con- 
tribution to the static shift S^lj stems from the diagonal 
terms = ly' in Eq. ([T3]), whereas the non-diagonal terms 
are suppressed by a factor uj'^tl^, , respectively. 

We henceforth work in the interaction picture with re- 
spect to the renormalized energies so that the density 
matrix element pjji reads 



P.y(i) = T7e-^-'(*^+^^^^'(*^. (16) 



The time-dependent phases (pjj' (t) — ipj (t) — ipj/ {t) decay 
to zero after a rather short time t ~ uj^^ and, thus, influ- 
ence the decoherence process only during a short initial 
stage. 

b. Damping factors. The coherence loss is given by 
the damping factors exjp[—Ajji(t)]. Inserting the ohmic 
spectral density (fTT|) into expression lU, we obtain for 
them the explicit form 



,{t) _ 



N 



iujctjj' 
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r'(fe^[i + -ct]) r^(fe^[i - ^^.t]) (1 + ^It-) 



r^fe^^Il + -c(i - )]) r2(^[l - ^u.^{t + 1,,)]) (i + (irllg^) 



(17) 



where F denotes the Euler Gamma function. The nomi- 
nator of the second factor in Eq. (fT7|) itself is already of 
physical interest: It describes the decoherence of a sin- 
gle qubit in the absence of the other qubits^— One can 
identify three stages in the single-qubit time evolution:^ 
Very shortly after the preparation, i.e. for times t < uj^^, 
the fluctuations of the bosonic field are not yet effective, 
leading to a "quiet regime" in which essentially no single- 
qubit decoherence takes place. At an intermediate stage, 
u}~^ < t < h/k^T^ the main origin of single-qubit deco- 
herence is vacuum quantum fluctuations. They lead to 
an initial slip of the coherence which we discuss in detail 
below. Finally for times larger than the thermal coher- 
ence time, t > h/k^T , thermal fluctuations dominate the 
coherence loss. The dephasing will finally be complete, 
i.e. a single qubit that starts in a superposition will loose 
all quantum coherence due to dephasing caused by the 
one-dimensional bathj^ 



For the spatially separated qubits prepared in the W 
state (fT2)) that we focus on here, the transit times tjji 
introduce additional time scales after which the denom- 
inator of the second factor in Eq. ([T7]) becomes rele- 
vant. Since the denominator ultimately decays equally 
fast as the nominator, decoherence will come to a stand- 
still. This interesting behavior does not occur for a single 
qubit and will now first be investigated for the case of a 
qubit pair. 



B. Entanglement of qubit pairs 

It is instructive to first consider the case of two qubits, 
N — 2, which represents the most basic system in which 
a spatial qubit separation infiuences quantum coherence. 
The robust initial state IT^ then reduces to the maxi- 
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mally entangled Bell state 

|W^2) = 4^(110) + 101)) = -i=(|l) + |2)), (18) 



V2' 



where the last expression refers to the basis (|13p . 

For a bipartite system, the degree of entanglement can 
be measured by the concurrence^ 



C = max 



{0, 



(19) 



where the Xi denote the eigenvalues of the matrix 
p(jyiUy2P*<Jyi<Jy2 in decreasing order and p* is the com- 
plex conjugate of p. For qubit pairs that are initially 
prepared in state ((T8)) and subject to pure phase noise, 
the concurrence can be expressed at all times by the 
absolute value of a single non-diagonal element, namely 
C{t) = 2|/9oi,io(i)| = 2|/Oi2(i)|, irrespective of the spatial 
separation. 

Figure [2] shows the time dependence of the concurrence 
of a qubit pair prepared in the robust entangled state 
(|18p for various spatial separations. For vanishing sepa- 
ration, we find the concurrence C(t) = 1, which means 
that the entanglement indeed remains perfect and justi- 
fies the designation "robust entangled state" ,2^21 This is 
tantamount to saying that the state IW2) is an element 
of a so-called decoherence-free subspace of the two-qubit 
Hilbert space. However, Fig. [2] also shows that the con- 
currence of the robust state does decay if the qubit sep- 
aration X12 is finite. The decay lasts until the transit 
time ti2 is reached. From then on, the concurrence re- 
mains constant, so that both the entanglement and the 
coherence become stable. Thus, we can conclude that 
for spatially separated qubits the state IVF2) is not an 
element of a decoherence-free subspace, but rather of a 
decoherence-poor subspace. Its emergence from the ex- 
act solution (fT7|) will be discussed for the more general 
case of N qubits in the subsequent section. Moreover, a 
more intuitive picture will be drawn in the framework of 
a causal master equation approach in Sec. IIVI 

Further information is provided by the value to which 
the concurrence saturates. Figure [5] shows that C{t 
00) is influenced by both the transit time between the 
qubits and the properties of the bosonic field. The exact 
value can be obtained from Eq. pT|) and reads 



( fcBT \ 
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(20) 



Figure [3] shows this final concurrence C(oo) as a function 
of the transit time ti2 for various temperatures T. As 
for the time-evolution, three regimes can be identified: 
For the (unphysically small) separations X12 < c/lUc, the 
concurrence remains at C — 1, while for c/coc < X12 < 
hc/k^T , the entanglement is no longer perfect, but still 
at an appreciably large value. For large separations, 
X > hc/k^T, the concurrence essentially decays to zero. 




FIG. 2: (Color online) Time evolution of the concurrence C 
for two qubits initially prepared in the robust state (|18|l for 
various transit times t\2'- exact time evolution (solid lines) 
compared to the results obtained from the causal master 
equation derived in Sec. IIVI (dashed). The temperature is 
fer = 10~*fiwc and the coupling strength a = 0.005. For 
Q, = lQ~^LUc, the time range corresponds to 40 coherent oscil- 
lations. Inset: Blow-up on a logarithmic scale for the transit 
time trz ~ 5 X lO^'/ojc. 
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FIG. 3: (Color online) Final value of concurrence for the 
robust state IW2) as a function of the spatial separation 
X12 ~ cti2 for various temperatures. The coupling strength 
is Q = 0.005. 



The latter limit is a prerequisite for the application of 
quantum error-correction schemes that assume that the 
qubits experience uncorrelated noise. 

Summarizing our two-qubit results, we find that the 
two-qubit concurrence can saturate at a stable value as 
long as the separation of the qubits stays finite, i.e. as 
long as a;i2 < 00. Remarkably, the coherence of a single 
qubit in the same environment (i.e. when only one qubit 
is present) decays completely. The entanglement dynam- 
ics that we find for the robust-entangled Bell state ^8]) 
is highly non-Markovian and it is obvious that such a 
behavior can not be described by a Markovian master 
equation such as the standard Bloch-Redfield equations. 
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as we discuss below. This fact, and in particular the final 
saturation of concurrence as shown in Figs. [5] and [3] con- 
tradict expectations formulated in Ref. uEh calculations 
in that work involve non-entangled initial states, while 
we study robust-entangled Bell states. 

Finally, we want to note that an entanglement satu- 
ration under pure dephasing may appear for baths with 
super-ohmic spectral densities as well )^^'^^i^° even in the 
limit of large qubit separations. This behavior is in con- 
trast to the intuition emerging from recent work,— 
namely that the entanglement would always decay to zero 
when qubits are coupled to independent heat baths. For 
super-ohmic baths in higher dimensions, a finite final en- 
tanglement does not necessarily stem from spatial cor- 
relations. This becomes obvious from the fact that for 
super-ohmic baths as for example a three-dimensional 
phonon field, the coherence of single qubits exhibits a 
similar saturation. 



C. iV-qubit fidelity 

Let us now turn to the intriguing question how the 
previous results can be generalized to arrays of qubits. As 
already described above, the focus will be on the scaling 
of the decoherence as a function of the system size N 
of linearly arranged qubits. In particular, we study how 
the amount of entanglement evolves for the qubits that 
start in the A^-qubit W state (|12p and, due to dephasing, 
at later times must be described by the iV-qubit density 
matrix 

With the exact dynamics known, the only remaining 
question is how to quantify the entanglement. If there 
were no interaction with the environment, then the qubits 
would remain in their pure entangled W state p2p . which 
in density-matrix notation reads p(0) = |WAr)(W^Ar|- Af- 
ter the dissipative time evolution, the qubit state deviates 
from this "ideal" output state pidced{t) = p(0). The ques- 
tion is how much. A proper measure for this quantity is 
the fidelity^ F(t) = tr{p(t)pidoai(t)}, which m our case 
reads 



F{t) =tr{piO)p(t)} = {WN\m\WN) 



(21) 



In general, the fidelity is bounded by < F < 1, where 
F — 1 corresponds to a pure stateM For qubits subject 
to pure dephasing, the somewhat more strict condition 
J2j P^ji^) < < 1 apphes, because the populations do 
not change. In particular the inequalities 1/A^ < F < 1 
will hold for the initial state jV^Af), as illustrated below. 
To give another argument why fidelity makes a good mea- 
sure of entanglement for our particular purpose, we em- 
phasize that for two qubits prepared in the state , the 
fidelity directly relates to the concurrence via the relation 
C{t) = 2F{t) - 1. 

We note in passing that for other initial states that are 
not pure A'^-qubit entangled states, one should be careful 
to use fidelity as an entanglement measure, for example 



because the fidelity may remain constant while the sys- 
tem undergoes nontrivial dynamics ^^-^ Finding other en- 
tanglement measures for three or more qubits is an active 
field of researchj22iSii^i^ Their numerical evaluation can 
be rather involved, especially for larger systems. These 
issues need not concern us here, since we start with an 
A^-qubit entangled pure state for which the fidelity is a 
good measure of entanglement. The fidelity has the ad- 
ditional advantage that it is easily evaluated analytically 
for larger systems as well. 

For the robust state p^ . the fidelity becomes 



1 ^ 



(22) 



where the coherences pjji are given in Eq. (|16p . The time 
evolution of the fidelity for Af = 3 and for A^ = 6 qubits is 
shown in Fig. |31 For low temperatures [panel (a)] , we find 
that the fidelity decay is slowed down whenever a transit 
time is reached, i.e. at times t = tjji . This resembles 
the behavior of the concurrence for two qubits shown 
in Fig. [51 For a larger number of qubits, the fidelity 
saturates at a lower value. 

In order to gain a more quantitative understanding 
of the fidelity saturation, we focus on the final fidelity 
F(oo). In the zero-temperature limit, which provides a 
lower bound for the coherence loss, we find from Eqs. 
and (|17p the density matrix elements pjjt (oo) = [1 + 
^^c^LO' - Hence the fidelity ^ saturates 

to 



F(oo) 



1 

TV 



N-l 



l-q/N 



(23) 



Although we start out with a "robust" entangled state, 
this final fidelity F{oo) can be as low as 1/A^, which 
marks the large-distance limit ti2 oo. More gener- 
ally, we find that for zero temperature, the final fidelity 
decreases with increasing cutoff frequency tUc , increasing 
spatial separation, and for a larger qubit-bath coupling 
strength a. 

An intriguing aspect of the fidelity is its scaling be- 
havior as a function of the system size A^. Will F{oo) 
decay to zero for larger arrays of qubits, or converge to 
a finite value? For large A^ we can neglect the term 1 /N 
in Eq. (|23p and replace the sum over q by an integration 
over the continuous variable x = q/N. Then we obtain 



Fi^) 



dx 



1-x 



(1+A^2^2i2^^2)4c 



22F1 ( i,4a,^,-ArVt' 



c ''12 
2, ,2-f2 \l-4a 



1 - (1 + N^UJ^JI^) 

7VVi?2(l-4a) 



(24) 



where the evaluation of the integral yields Gauss' hy- 
pergeometric function 2^1- The expression (I24p is valid 
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FIG. 4: (Color online) Exact time evolution of the fidelity 
F{t) (sofid lines) and the result obtained from the causal mas- 
ter equation (dashed lines) for A'^ = 3 and N = 6 qubits, re- 
spectively. The temperatures are ksT — 10~^ /huic (a) and 
ksT — 10~'^ /hujc (b). As in Fig. (2] the transit time between 
nearest-neighbor qubits is ti2 = 5 x W^/ujc and the qubit-field 
coupling strength a = 0.005. The inset in panel (b) shows the 
data for A'' = 3 on a logarithmic time axis. 



for general coupling constant a, but beyond weak cou- 
pling its value is rather small. We will now approximate 
Eq. for a ^ 1. Furthermore, usually many cut- 

off wavelengths 2t:c/lUc will fit between two neighboring 
qubits, so that Wcii2 ^ 1, as we argued above. Since we 
already assumed iV ^ 1 to arrive at the integral ([M)) . 
we are surely in the limit Nu!cti2 2> 1. Then we can ap- 
proximate the hypergeometric function by its asymptotic 
expansion for large fourth argument. We finally obtain 



F(oo) 



r(i-4a) 1 



r(f-4Q;) l-4a 



-8a 



(25) 



Hereby we found the important result that although the 
final fidelity is smaller for larger systems, the scaling is 
only algebraic in N . This stability under dephasing is a 
property of the initial AZ-qubit W state (fT2|) . Clearly, 
for nonvanising ti2 this state lives in a decoherence- 
poor rather than in a decoherence-free subspace. Equa- 
tion ((25)) shows that at zero temperature, the final fi- 
delity is determined by two dimensionless numbers, the 
one number being a and the other the ratio between the 
array length Ncti2 and the cutoff wavelength 27rc/wc. 
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FIG. 5: (Color online) Final fidelity (|23|l as a function of 
the number A'' of qubits at zero temperature for two cou- 
phng strengths a = 0.001 (blue circles) and a = 0.01 
(red diamonds). The nearest-neighbor transit time is again 
<-L>cti2 = 5 X 10*. The crosses (orange) and plus signs (green) 
mark the approximative result (|25p , respectively. The dashed 
line indicates the lower bound 1/JV of the fidelity. 



In Fig. O we compare for two values of a the exact 
expression ([25)1 for the final fidelity as a function of N 
with the weak-coupling approximate result (j25p . Clearly, 
for state-of-the-art well-isolated qubits with typically a — 

0. 001, the agreement is excellent already for iV > 5, while 
for a = 0.01 convergence is reached for > 10. The 
figure clearly shows that in the weak coupling limit a <IS; 

1, the final fidelity F{oo) in (PS)) is almost independent 
of the length of the array. For a < 0.005 the factor in 
square brackets in Eq. ([^5]) is less than 1.05, so that the 
large- A^ expression for the final fidelity could be further 
simphfied as F{oo) ~ {Ncjcti2)~^" ■ 

The above estimates were derived in the limit of 
strictly zero temperature, so that the question arises up 
to which temperature they still represent a reasonably 
good approximation. A closer inspection of the exact 
result (fT7)) reveals that this is certainly the case if the 
condition 



he 



> Nch 



(26) 



holds, i.e. if the thermal coherence length of the bath 
is much larger than the length of the array. Assuming 
T = 10 mK, the thermal coherence length is 2.3 /im. This 
value corresponds to an array length of 24 qubits with a 
nearest-neighbor distance Xi2 = cti2 = 100 nm. The final 
fidelity that we obtain for T = in Eq. (|25p can therefore 
be considered as an upper bound for what could be real- 
ized in state-of-the-art quantum information processing 
experiments on arrays of qubits. 



IV. CAUSAL MASTER EQUATION 



In Sec. Illli we found that the analytical solution for 
the qubit dynamics can involve rather complex expres- 
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sions and, thus, an intuitive picture of the observed be- 
havior can be hard to find, even though the exact solu- 
tion is known. Thus for a more quahtative understand- 
ing, one can benefit from an approximative treatment in 
the spirit of a Bloch-Rcdficld master equation approach. 
Moreover, such an approach enables a symmetry analy- 
sis of the dissipative time evolution. This can provide 
additional insight in cases in which tracing out the bath 
degrees of freedom reveals symmetries that are obeyed 
by the dissipative central system, but not by the system- 
bath Hamiltonian. 

Bloch-Redfield equations are based on a perturbative 
treatment of the qubit-environment coupling, followed 
by neglecting memory effects in the kernel of the re- 
sulting quantum master equation. Thereby one entirely 
ignores the dependence of the dynamical equations on 
the qubits' history and, thus, on the initial preparation. 
The resulting master equation then assumes the struc- 
ture dpjn,n/dt = J2in',n' ^mnm'n' Pm' ,n' ■ Such CquatioUS, 

however, fail to reproduce the stepwise decay of the con- 
currence observed in Sec. lIIII For the concurrence C{t) of 
the robust Bell state (|18p . they even yield dC/dt = for 
all distances X12, in clear contrast to the exact solution. 
In other words, strictly Markovian master equations can 
predict spurious decoherence-free subspaces. In this sec- 
tion, we shall derive a master equation that does not 
suffer from such a shortcoming. 



A. Markov approximation and beyond 

Taking the above considerations as a motivation, we 
now derive a generalization of the Bloch-Redfield mas- 
ter equation that is able to capture the retardation ef- 
fects stemming from a finite sound velocity. In doing so 
we pursue closely the standard approach that leads to 
a Markovian equation of motion. In this way, it can be 
identified where it fails and how to improve it accord- 
ingly. 

Starting with the Liouvillc-von Neumann equation ([7]), 
we employ a projection-operator formalism to formally 
eliminate the bath degrees of freedom. In second-order 
perturbation in the qubit-bath coupling, we obtain for 
the reduced qubit density matrix the master equation 



dt 



N 

E 

um' —1 



dr Ky^i{t, t) p(t) 



where t > to and the super-operator 
K..'(i,T)[-] =5...(r)[X,(i), [X,,{t-T), •]] 



(27) 



(28) 



with the anti-commutator {•,•}. For a sketch of the 
derivation, see Appendix [Bj Although this master equa- 
tion is in a timc-convolutionless form, it still is non- 
Markovian owing to the explicit dependence on the initial 
time tf). The integral kernel features the real part 



(t) and the imaginary part A^,^' (r) of the bath corre- 
lation function t'rh[^u{T)^u' P}^]- For our bath model they 
can be evaluated explicitly and read 

Su.' (r) = \[SiT- U,, ) + 5(t + U,, )] , (29) 

A,.' - ) + -A(t + t,,, )] , (30) 



where 



^{'^) = / dio J{lo) cos(ti)i) coth ( — 

Jo X^kbT 
1*00 

A{t) = dio J{uj) sm{ujt) , 

"'0 



(31) 
(32) 



are the usual symmetric and anti-symmetric bath corre- 
lation functions respectively. 

We first consider the local terms, i.e. those with 1/ = v', 
for which the time shift 1,^1,1 vanishes. Then the correla- 
tion functions ^ and ^ reduce to Eqs. ^ and ([5^ . 
respectively, and one can introduce a Markov approxima- 
tion in the usual way: If the correlation functions S{t) 
and A{t) contribute to the integral in Eq. essentially 
in a small time interval of size Tb around t = 0, then for 
t — to ^ Tb , we can extend the r-integration to infinity, 
i.e. we set 



dr Huuit,T) 



dr Ky^{t,T) . 



(33) 



This expression implies a coarse-graining in time so that 
the resulting master equation is valid only for time steps 
not smaller than the bath correlation time Tb. In general, 
the bath correlation time depends on the properties of the 
spectral density J {to) and the temperature of the bath. If 
the temperature is not too low and the spectral density 
is fairly smooth and decays sufficiently fast for lu —^ 
and w — !■ 00, as is the case here, then the correlation 
time is only weakly temperature dependent and reads 

Tb « iK-^ 

In the above treatment of the local terms — ly' , we 
have followed the route towards a Markovian equation of 
motion. For the nonlocal terms, however, the arguments 
of the last paragraph are no longer valid. For v ^ v' 
the correlation functions ((29|) and pO|) are not peaked 
at T = 0, but at T = ij,^/ which for a realistic qubit 
separation typically exceeds the bath correlation time Tb. 
Then we have to distinguish the cases t — to < Tb and 
t— io > Tb- In the former case, the peak of the correlation 
functions lies outside the integration interval so that the 
integral is small and, consequently, will be neglected. In 
the latter case, by contrast, the peak fully contributes so 
that the integral can again be extended to infinity. In 
summary, this means 



t-ta 



dTK^l,\t,T) 



10 



(34) 
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where Q{t) is the Heaviside step function. For ly = i/', 
this expression coincides with Eq. ([55]) . 

Inserting the approximation ([34|) into the weak- 
couphng master equation (|27p and setting again the ini- 
tial time to = 0, we find the causal master equation 
(CME) 



dt 



Pit) = i?(t)/5(i) 



(35) 



with the time-dependent superoperator 



R{t)[-] = 



(36) 



The step functions ensure causality which requires that 
the cross terms can only be active after the propagation 
time between the respective qubits has passed. In the 
limit of vanishing separation, the causal master equation 
reduces to a standard Bloch-Redfield equation. In the 
following, we will demonstrate that the causal master 
equation reproduces the results of Sec. |llT] rather well, 
while a standard Bloch-Redfield approach clearly fails. 



B. Master equation for pure dephasing 

Let us now apply the causal master equation to 
the problem defined in Sec. |TT] and test the results 
against the exact solutions presented in Sec. IIIII Since 
— commutes with the free Hamiltonian Hq^ 
the interaction-picture coupling operators stay time- 
independent, Xi,(t) = a^z- Then the time integration in 
the causal Bloch-Rcdfield tensor involves only the 
bath correlation functions and wc obtain for the density 
matrix element Pm.n the equations of motion 



d 
dt 



-A™^(0 + *&"(i) 



CME/ 



/^m,n 



(37) 



with the damping 



N 



1^,1/'— 1 

X [(-1)""+'".^ + (-1)"-+". 

and the phase shift 

N 



(38) 



iCME 



aujc 



1 

|^(^_X)"'"+™'^' — (— 1)" 



— 1 

X 



(39) 



This master equation is non-Markovian due to the ap- 
pearance of the step functions, which change the effec- 
tive damping and phase shift whenever a transit time 
t^^i is reached. These stepwise time-dependent frequency 
shifts and decay rates are characteristic features of the 
causal master equation. As we will see, because of these 
steps the CME follows more closely the smooth time- 
dependent variations of shifts and decay rates of the exact 
dynamics than the standard master-equation formalism 
manages to do with its static shifts and decay rates. 

As in Sec. IIIII we first consider as an initial state the 
robust Bell state (|18p . From the master equation (|37p . 
we find that the concurrence Ccme = 2|poi,io| obeys 



dt 



C. 



CME 



[i-e(t~ti2)] c\ 



CME ■ 



(40) 



This differential equation is readily integrated to provide 
the solution 



Ccme(0 



(41) 

i.e. the concurrence decays exponentially until the tran- 
sition time is reached and thereafter remains constant. 
This clear separation of two dynamical regimes facilitates 
an intuitive interpretation: For times t < tu, the qubits 
have not "seen" each other, so that we are in a regime of 
single-qubit decoherence. Indeed, during this first time 
interval the causal master equation ([37]) coincides with 
a standard Bloch-Redfield approach in which the qubits 
are coupled to independent heat baths. Consequently, 
the relative phase between the qubits is randomized and 
the concurrence decays. However, for t > ti2, both qubits 
experience correlated quantum noise and undergo collec- 
tive decoherence. Thus, the concurrence decay comes to 
a standstill and a decoherence-poor subspace can emerge. 

This time evolution is compared to the exact solutions 
in Fig. [21 We find that generally the causal master equa- 
tion describes the slow decay of the concurrence and its 
saturation very well. At very short times, however, the 
causal master equation does not capture the initial slip of 
the concurrence. The reason for this is, that the dynam- 
ics on timescales that are comparable to the bath correla- 
tion time cannot be resolved in a coarse-grained time ap- 
proximation underlying the causal master equation. The 
same generally holds true for Markovian quantum master 
equations, in particular for the standard Bloch-Redfield 
equation. At larger times, the benefits of the causal mas- 
ter equation become obvious: Since the Bloch-Redfield 
treatment is recovered by setting the transit time ti2 — 0, 
Eq. (|40p reveals that dCBR/dt — for all times, i.e. the 
concurrence remains at its initial value. This spurious ro- 
bustness of the concurrence is, however, in clear contrast 
to the exact result. For pure phase noise, the popula- 
tions of the system eigenstates are conserved and, conse- 
quently, the final state is generally unrelated to thermal 
equilibrium. Therefore, the question whether the master 
equation can describe the grand canonical ensemble of 
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FIG. 6; (Color online) Difference between the exact result 
and the causal master equation result (|4ip for the final value 
of the concurrence, 5C = C(oo) — Ccme(oo), for the robust 
entangled 2-qubit state (|18p . The coupling strength is a = 
0.005. The inset depicts SC for the fixed temperatures 9 = 
kBT/huc = lO""* (solid), e = 10"^ (dashed) and = 10"^ 
(dash-dotted), respectively. 



the system coupled to the bath is in the present context 
ill-posed. 

For a quantitative investigation of the quahty of the 
causal master equation, we compare the final values of 
the concurrence, limt^oo C{t). Figure [6] depicts the dif- 
ference of the exact solution and the causal master equa- 
tion result, SC — C(oo) — Ccme(oo), as a function of 
transit time and temperature. For realistic temperatures 
T < 10~^hujc/kB, Fig. [6] shows that the final value ob- 
tained from the causal master equation exceeds the exact 
result, so that Ccme(oo) provides an upper bound for the 
concurrence. In particular, for lOOH/kBT < ti2 < i^^^^ 
the agreement is almost perfect. Note that for very high 
temperatures T > lO^^ftwc/^B, the difference SC as- 
sumes also positive values. 

Let us finally apply the master equation ((37)) also to the 
case of a linear A'^-qubit arrangement with equal nearest- 
neighbor spacings Xi2 as discussed in Sec. lIII Cl We again 
consider the initial preparation in the robust state (jl2p 
using the shorthand notation (fTS]) . To calculate the fi- 
delity F defined in Eq. ([22]) . we need to compute the 
values of the density matrix elements pjj' . Unlike for the 
two-qubit concurrence, phase shifts now also play a role. 
The time-dependent phase shift (t)™^it) [see Eq. ([39]) ] 
can be written as the difference of two terms, the one 
only depending on j and the other only on j' . Both terms 
describe stepwise time-dependent frequency shifts of the 
corresponding qubits. Interestingly, after the longest 
transit time t > tiN, these shifts in the CME become 
static and (j)'^Jf^{t) agrees with the exact result P^ . i.e. 
(f>^p^^(t) = Silj — Srijt . For an unambiguous comparison 
of fidelities in the exact and in the causal master-equation 



formalism, it is an important result that we can work in 
the same interaction picture with the same renormalized 
frequencies flj Qj + Sflj . It remains to be discussed 
what is the effect of the stepwise frequency shifts in the 
CME for times t < ti^. It is the non-diagonal terms 
V ^ v' in Eq. ((39|) that make up the difference between 
the frequency shifts at time t — Q and the exact static 
renormalization at times t > tiN. However, this differ- 
ence is very small due to large factors w^i^^, 3> 1 in the 
denominator of Eq. (|39p . and can safely be neglected in 
the following. 

From the causal master equation (|37p we then obtain 



d ^ / \ 



[l-e(i-i,y)]p,y(t). (42) 



Interestingly enough, in the present case only two types 
of terms of the master equation ([37]) contribute: the local 
terms with v ~ u' and those with |i/ — = |j — As 
a consequence, the decay rate of pjj' (t) changes only at 
the transit time tjj' . 

In order to evaluate the fidelity ([22ll . we integrate 
Eq. (|42p and sum over all density matrix elements pjji {t) , 
so that we obtain 



N 



j,3' = l ^^'^^ 



For vanishing qubit separation, the master equation pre- 
dicts FcME{t) = 1, i.e. a decoherence-free behavior. For 
Xi2 > 0, however, all coherences pjji initially decay and 
so does the fidelity. When the smallest transit time ti2 
is reached, a perfect correlation between nearest neigh- 
bors is built up and the coherences Pjj+i saturate. Since 
these iV — 1 coherences are no longer time-dependent, the 
fidelity decay is reduced accordingly. This process con- 
tinues until ultimately all transition times have passed, 
i.e. until t = fijv, and the fidelity decay comes to a stand- 
still. 

In Fig. m we compare this behavior to the exact so- 
lution. For very low temperatures [panel (a)], we ob- 
serve that the master equation reproduces the reduction 
of the fidelity decay whenever a transit time is reached. 
The relative difference between the exact result and the 
causal master equation is of the order of 10% as in the 
case of two qubits. If the temperature becomes larger, 
so that the nearest neighbor separation exceeds the ther- 
mal coherence length [see Fig. [3] and the discussion after 
Eq. (|20p]. then the fidelity decay is determined by ther- 
mal noise. In that case, F already saturates to a rather 
small value before the first transit time is reached. The 
inset of Fig. shows that at very short times when 
the decay sets in, i.e. when the CME is not yet differ- 
ent from a standard master equation for independently 
dephasing qubits, that then the CME result can devi- 
ate significantly from the exact result. Again we wish 
to stress that this deviation is due to the coarse-grained 
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time approximation that is inherent in a standard master 
equation approach as well. However, in contrast to the 
latter, the CME and the exact results agree very well at 
later times. 



V. DISCUSSION AND CONCLUSIONS 

Decoherence of an array of qubits can be a rather com- 
plex process that proceeds in several qualitative different 
stages, even in the case of pure phase noise which we in- 
vestigated here: If the qubits are coupled to the bosonic 
field of the substrate vibrations, the dynamics during a 
first, very short period is essentially noiseless. In a second 
stage, the bosonic vacuum fluctuations arc most relevant, 
while finally, thermal noise dominates. The spatial sepa- 
ration of the qubits brings in further time scales, namely 
the transit times of sound waves between the qubits. 

Here we considered robust entangled A^-qubit W 
states, which do not decohere for vanishing qubit sep- 
arations. For finite separations we find instead that 
the qubits start to dephase. However, the dephasing 
slows down whenever the elapsed time reaches a tran- 
sit time, until it eventually comes to a standstill. The 
final A'^-qubit quantum coherence increases with decreas- 
ing qubit-qubit separation, qubit-bath coupling strength, 
cutoff frequency, and temperature. By contrast, single 
qubits in the same one-dimensional environment would 
loose all quantum coherence for all finite values of these 
bath parameters. Note that, the two-qubit W state is 
identical to the robust-entangled Bell state. The satu- 
ration of entanglement that we find does not occur for 
fragile two-qubit statesi^iS? 

Cooperative effects can be advantageous or detrimen- 
tal, depending on the specific protocol that one has in 
mind. For example, one may fight decoherence by creat- 
ing decoherence-free subspaces. To that end, one could 
bring the qubits close together and use the W states for 
quantum information processing because their entangle- 
ment is robust. Nevertheless, the qubits must be suffi- 
ciently well separated, either to enable their individual 
manipulation or because of their finite extensions. We 
found that this requirement prevents the realization of 
decoherence-free subspaces. Cooperative effects are still 
advantageous, since decoherence-poor subspaces are built 
up instead. Good results require the length of the array 
of qubits to be smaller than the thermal coherence length 
hc/k^T. 

Alternatively, one may wish to implement active quan- 
tum error-correction schemes, where logical qubits are re- 
dundantly encoded into several physical qubits. Here, by 
contrast, the cooperative effects are detrimental, since 
standard error-correction schemes^^iS and recent gen- 
eralizations to non-Markovian baths^i^ will only work 
perfectly if the physical qubits couple to spatially uncor- 
related baths. Our discussion demonstrates that neigh- 
boring physical qubits should then be separated by more 
than the thermal coherence length hc/k^T. By reducing 



the temperature, single-qubit dephasing is suppressed, 
but the assumption of uncorrelated baths becomes worse. 
This suggests that there may be an optimal working tem- 
perature for quantum error correction, given a geometry 
of physical qubits. 

Thus our calculations show how well decoherence-free 
subspaces or quantum error correction-protocols could 
be realized with linear arrays of qubits. The aforemen- 
tioned conflicting requirements for both strategies seem 
to rule out the implementation of both strategies in one 
experiment. 

One should keep in mind that bit-flip noise, which was 
not considered here, will reduce coherence further, al- 
though typically on a longer time scale. Including bit- 
flip noise usually renders the qubits-environmcnt model 
no longer exactly solvable. Therefore one has to resort 
to approximation schemes like, e.g, a master-equation 
approach. As we have shown, the common Bloch- 
Redfield master equation cannot account for the intrin- 
sically non-Markovian effects stemming from the spatial 
separations. For the present model, the Bloch-Redfield 
approach predicts spurious decoherence-free subspaces, 
which in fact are at best decoherence-poor. Very re- 
cently, it was found^ that for bit-flip noise due to a 
homogeneous isotropic Markovian and i/iree-dimensional 
reservoir, there is no multi-particle decoherence-free sub- 
space outside the Dicke limit, i.e. whenever the qubits are 
not co-located. The situation may be different for one- 
dimensional structures for which we find that at least for 
pure dephasing, decoherence-free subspaces can become 
imperfect if the qubits are separated. 

In order to capture delocalization effects with a mas- 
ter equation, we have derived a modified Bloch-Redfield 
approach that ensures causality for the qubit-qubit in- 
teraction mediated by the substrate. It proved to be 
reliable for parameters for which the standard master 
equation for a single qubit is reliable. This is the case 
for sufficiently high temperatures or small enough cou- 
pling strengths such that initial-slip effects are small. 
A characteristic feature of the proposed causal master 
equation is that it selects the Bloch-Redfield kernel de- 
pending on the time elapsed since the preparation. This 
means that the time evolution is governed by a time- 
dependent Liouville operator which renders the dynamics 
non-Markovian. Besides being a proper tool for studying 
retardation effects in models that do not possess an exact 
solution, the causal master equation describes the time 
evolution in an intuitive and concise manner. Thereby, 
it enables decoherence studies with algebraic methods 
which possibly will provide suggestions for coherence sta- 
bilization. 

In this sense, our studies of the interplay between pure 
dephasing and a spatial qubit separation can only be a 
first step towards a deeper understanding of such phe- 
nomena. In particular the inclusion of other quantum 
noise sources, which are also present in real experiments, 
will complement the picture drawn above. 
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APPENDIX A: EXACT REDUCED DYNAMICS 

In this appendix, we outline the derivation of the exact 
solution ([5]) for the reduced qubit dynamics for the case of 
pure dephasingj^i2ii^iS2i. In the usual interaction picture 
with respect to the uncoupled qubits and the bath, the 
coupling Hamiltonian ([3]) reads 

H^y,it)^V{t) + V\t), (Al) 
with the interaction 

V{t) = n5]a,,^gfc6fce^('=-°--'=*) . (A2) 

k 

Our aim is to calculate the time evolution of the reduced 
density matrix p{t) = trb R{t) = tr\^[U{t)R{Q)W {t)] gen- 
erated by the propagator 



U{t)=T exp ( - / dsH^^{s) 



(A3) 



where T denotes the time-ordering operator. In a first 
step, we find [V{t),V'^ {t')] =J{t-t'), where 



(A4) 



k w' 



Since ^{t),Y(t')\^^\t),V\t')\=Q and \j{t),V{t')\ 
= [/(<), V"^ (/:')]= for all times t and t', we can use 
the Baker-Campbell-Hausdorff formulaic to express the 
time-ordered exponential (|A3|) as 

U(t) = exp 1 4 / dsHqi,{s) f ds f ds' 

I «^ Jo Jo Jo ^^5^ 

xfis-s') [e{s - s') - 0is' - s)]' 
The first term in the exponent can be written as 

expj^^ 7?qb(t)ds| =PJi?fc(^(7,,y,fc) , (A6) 

where we defined y^k = gke~'^''^"{l — e*'^'"*)/a;fc and the 
displacement operators Dk{X) — exp{Xbl,~ X^'bk}- The 
second term in the exponent is a c-number and provides 
the time-dependent phase factor exp[i(/)(t)] with 

^2^ 2^ 9k ^ -(JuzCT,,,e''''=-"' . (A7) 



So far we have found for the propagator the expression 
U{t)^X{ Du (T..y.fe) e^*(*) . (A8) 



For the factorizing initial condition i?(0) — p{0)p^, the 
matrix elements of the reduced density operator become 

Pm.„(t) = trb(m|C/(t)p(0)p;;'iC/t(i)|n) 

= Pm,n(0) trb {pI^ (n| C/t (t) |n) (m|C/(i) |m) } . 

(A9) 

In the second equality of Eq. (|A9|) we have used the cyclic 
property of the trace and the fact that the computational 
basis elements |n) are eigenstates of the propagator. In- 
serting the propagator (jASP into (jA9p and assuming an 
isotropic coupling g^k — 9k, we find for density matrix 
element /5m,n(^) the phase 



0m,nW = (m|0(i)|m) - (n|0(i)|n) 



(AlO) 



N 



EE 

k u^y' 



^ujkt ~ sm{ujkt) 

9 k —2 



(All) 



which in the continuum limit becomes the phase (jlOp . 
For the calculation of the remaining contributions in 
(|A9[) . we employ the relations Dj.{X) = Dk(—X) and 

Dk{X)Dk{Y) = e(^^'-^'^)/2^fc(A + Y), (A12) 

which hold for any commuting operators X and Y . Then 
we obtain 

Pm,n(i) 

= Pm,n(0)e^^-"W (A13) 
xtrb{pr(nin^UE^-y-0|n) 



(m| Yl Dk cFy^y^k) |m) | 



= Pm.n(O) e 



„(*)+>)„,„(*)] 



(A14) 



n trb {pllDk ( [(-1)"^ - (~1)"1 y,.k) } . 



An additional phase ?/m,n(0 stems from the commutator 
of the displacement operators Dk in (|A13p [see Eq. ljA12p ] 
and reads 



N 



??m,n(t) = ^ gl 



2 1 — COs(wfct) 



k uy=i 



(A15) 

X (-1)".^'+™^ sin(fcw) 

and vanishes for isotropic coupling g^k — 9k which we 
assume herein. 
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Finally, we have to evaluate the trace in Eq. (|A14[) . 
This can be accomplished conveniently in the basis of 
the coherent states \(3k), in which the equilibrium bath 
density operator © reads 



cq 
Ph,k 



1 



Tr{nk) 



(A16) 



The integration is over the whole complex plane and 
{nk) = [exp(ftti;fc/fcBr) — 1]^^ is the Bose distribution 
function. Pk are the complex eigenvalues of the annihi- 
lation operator hk^ i.e. bk\l3k) = Pk\Pk)- After inserting 
expression (IA16[) into Eq. (|A14[) . we integrate for each 
mode k over the complex plane. We finally end up with 



/5m,n(i) = /5m,n(0)e 



Z0n 



.(t)-A„.„(t) 



(A17) 



where 

Am,„(t) 



Hal 

k 



1 — COs(wfct) 



coth 



f^k 
2kBT 



N 



[(-1)™- - (-l)"'']e*'=^^ 



(A18) 



Note that Am,n(0 is real- valued and thus accounts for 
the damping of the matrix element. In the continuum 
limit, it assumes the form ©. 



{R{t)) of the density matrix, which describes the reduced 
dynamics. The projection operator (|B2p is constructed 
such that the factorizing initial condition ([5]) provides 
the identity {R{to)) = R{to). This allows one to write 
the formal solution of the dissipative system dynamics in 
the form 

(Rit)) ^{U{t,to))R{to). (B3) 
From Eq. (|B3[) . one can derive a formally exact 
quantum master equation for {R{t)) in two ways: 
The first possibility is the Nakajima-Zwanzig projector 
technique^i^ which leads to an integro-differential equa- 
tion for (R) which is nonlocal in time. The second possi- 
bility, on which we will focus here, is termed "time con- 
volutionless projection operator technique"— i^i^ and 
leads to a equation of motion of the form 



dt 



m)) = mm)) , 



(B4) 



which is local in time, but is governed by an explic- 
itly time-dependent super-operator IC{t). Although this 
equation possesses an apparently simple form, it gener- 
ally cannot be solved exactly and, thus, one has to resort 
to a perturbative treatment. In doing so, we expand 
the generator IC(t) in powers of the qubit-bath coupling 
parameter A, i.e. /C(t) = X]ri'^"^"(^)- Then a formal 
integration of Eq. (jB4p results in 



APPENDIX B: WEAK-COUPLING MASTER 
EQUATION 

For notational convenience, we write the Liouville-von 
Neumann equation ([7]) with the super-operator • ] = 
X[Hq\y{t), ■ ]/ifi, so that it reads dtR{t) — C{t)R{t), where 
A will serve as an expansion parameter. The formal so- 
lution can be written as R{t) = U{t,to)R{to), where 



U{t,to) = 



C(ti)dti 



u 
t 

to J to 



C{ti)C{t2)dtidt2 + 



(Bl) 



Moreover, we introduce the projection operator 
(•>-trb{-}®Pr 



(B2) 



which projects an operator of the full Hilbert space (of 
the qubits and the bath) to a direct product of a system 
operator and a time-independent reference state p'^ of 
the bath. Here, we are interested in the relevant part 



m)) = 



t / rt 

2 I 



1 + A / ICi{ti)dti + yl / K2{ti)dti 

to V-'to 



t ptl 



K.i{h)K.i{t2) dt2dti + 



R{to). 
(B5) 



to to 



Comparing the powers of A of Eqs. ()B3P and (|B5p . we 
find the relations 



/Ci(i)-(£(t)) 



(B6) 



/C2(i)= / {C{t)C{h)) ^{C{t)){Cih))dh, (B7) 

Jto 

so that the master equation reads 



-{R{t)) = {C{t)) + dh{C{t)C{h)){R{t)) . 

(B8) 

By inserting the definitions of the projection operator 
(jH2| and the LiouviUian C, we obtain Eq. ([27|. Note 
that the coupling ^ results in {C{t)) = 0. 
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